Hybrid method for understanding black-hole mergers: Head-on case 
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Black-hole-binary coalescence is often divided into three stages: inspiral, merger and ringdown. 
The post-Newtonian (PN) approximation treats the inspiral phase, black-hole perturbation (BHP) 
theory describes the ringdown, and the nonlinear dynamics of spacetime characterize the merger. 
In this paper, we introduce a hybrid method that incorporates elements of PN and BHP theories, 
and we apply it to the head-on collision of black holes with transverse, anti-parallel spins. We 
compare our approximation technique with a full numerical-relativity simulation, and we find good 
agreement between the gravitational waveforms and the radiated energy and momentum. Our 
results suggest that PN and BHP theories may suffice to explain the main features of outgoing 
gravitational radiation for head-on mergers. This would further imply that linear perturbations to 
exact black-hole solutions can capture the nonlinear aspects of head-on binary-black-hole mergers 
accessible to observers far from the collision. 

PACS numbers: 04.25.Nx, 04.30.-w, 04.70.-s 



I. INTRODUCTION 



Even prior to the complete numerical-relativity simula- 
tions of black- hole-binary mergers [ll-Q : black- hole colli- 
sions were thought to take place in three stages: inspiral 
(or infall), merger, and ringdown. During inspiral, the 
speed of the holes is sufficiently low and the separation 
of the bodies is large enough that the system behaves 
like two separate particles that follow the post-Newtonian 
(PN) equations of motion. Eventually, the black holes be- 
come sufficiently close that the dynamics given by the PN 
expansion significantly differ from those of full relativity. 
This stage is the merger, during which the two black holes 
become a single, highly distorted, black hole. The merger 
phase is brief; the strong deformations lose their energy 
to gravitational radiation, and the system settles down 
to a weakly perturbed, single black hole. The ringdown 
phase describes these last small oscillations of the black 
hole. 

Of the three stages of binary-black-hole coalescence, 
merger remains the most inaccessible to analytical tools. 
Nevertheless, full numerical relativity is not the only 
technique to have success investigating merger. Histori- 
cally, most analytical investigations of the merger phase 
arise from trying to extend the validity of perturbative 
techniques, particularly black-hole-perturbation (BHP) 
and PN theories. Those researchers working from a BHP 
approach try to push the approximation to hold at ear- 
lier times, whereas those employing a PN method at- 
tempt to stretch the technique to hold later into merger. 
Alternatively, one can see if there are exact, nonlinear 
analytical models whose dynamics can represent various 
aspects of black-hole-binary mergers. RezzoUa, Macedo, 
and Jaramillo |5| recently took this latter approach in 
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their study of anti-kicks from black-hole mergers. In their 
work, they showed that they could relate the curvature 
anisotropy on the past apparent horizon of a Robinson- 
Trautman spacetime to the kick velocity (computed from 
the Bondi momentum). Through appropriate tuning of 
the initial data, they could recover kick velocities found in 
numerical-relativity simulations of unequal-mass, inspri- 
aling black holes. While this type of approach is interest- 
ing and has proved successful, our work focuses on using 
perturbative approaches (and we will, therefore, take a 
more comprehensive look at the prior use of perturbative 
methods to understand mergers). 

From the BHP side, Price and Pullin [Q], initially, and 
many collaborators, subsequently, (see, e.g., [31 ) devel- 
oped the "close-limit approximation" (CLA). This tech- 
nique begins with initial data containing two black holes 
that satisfy the vacuum Einstein equations and splits it 
into a piece representing the final, merged black hole and 
perturbations upon that black hole. The exact form of 
the initial data varies in the CLA, but for head-on colli- 
sions of black holes, it typically involves some variation 
of Misner @ , Lindquist Q , or Brill-Lindquist ^^ time- 
symmetric, analytic, wormholc-like solutions. To investi- 
gate the late stages of an inspiral, the CLA usually begins 
from non-timc-symmctric, but conformally flat, multiple 
black- hole initial data set forth by Bowen and York [ll[ . 
Independent of the initial data, however, the CLA trans- 
lates the original problem of merger into a calculation 
involving BHP theory. The CLA does not allow for a 
very large separation of the black holes; as a result, only 
the very end of the merger is captured in this process. 
Moreover, "junk radiation" appears in the waveform be- 
cause the initial data do not describe the binary black- 
hole-merger spacetime in both the wave zone and the 
near zone. Unlike in full numerical-relativity simulations 
where the junk radiation leaves the grid during the well 
understood inspiral phase, in the CLA the junk radiation 
appears during the merger stage. This radiation (both 
from the absence of waves in the initial data and from 



errors in the near-zone physics), therefore, is difhcuh to 
disentangle from the physical waveform. 

The Lazarus Project (see e.g. [l2|) followed roughly 
the same approach as the CLA, but it used even more 
realistic black-hole-binary initial data for its CLA calcu- 
lation; namely, its initial data came from a numerical- 
relativity simulation just prior to merger. At the same 
time, however, because the initial data is now numeri- 
cal, one loses the analytical understanding of the proper- 
ties of spacetime near merger. More recently, Sopuerta, 
Yunes, and Laguna [l3| applied the CLA in combina- 
tion with PN flux formulae to obtain an estimate of the 
gravitational recoil from unequal-mass binaries (includ- 
ing binaries with small eccentricity [ij] ) . They proposed 
using more realistic initial data in the CLA, which Le 
Tiec and Blanchet [1^ ultimately carried out. Le Tiec 
and Blanchet chose to use the 2PN metric (keeping only 
the first post-Minkowski terms) as initial data for the 
CLA, and they applied it with considerable success to 
inspirals of unequal-mass black-hole binaries in a paper 
with Will [13 • Despite the improved initial data, this 
approach does not eliminate the problem of junk radi- 
ation discussed above. It would be of interest to see if 
even more realistic initial data, such as that of Johnson- 



McDanicl, Yunes, Tichy, and Owen [17[ would lead to 
improved results within the CLA. 

From the PN side, Buonanno et al., [I^, as well as 
Damour and Nagar [1^ take a different approach to un- 
derstanding the physics of merger. Using the Effective- 
One-Body (EOB) method [201, ^^^Y study the dynamics 
of the system until roughly the beginning of the merger 
phase. To obtain a complete waveform, they attach 
a ringdown waveform by smoothly fitting quasinormal 
modes to the EOB inspiral and plunge waveform. When 
they calibrate the two free parameters of this model to 
numerical-relativity waveforms, the EOB results match 
numerical-relativity waveforms precisely. In this method, 
they fit the PN dynamics and the ringdown at the light 
ring of the EOB particle motion; it is not immediately 
apparent, however, what this feature tells about the na- 
ture of merger. 

As a result, there remains a need to develop simple ana- 
lytical models that help reveal the behavior of spacetime 
during merger. Toward this end, it is helpful to delve 
deeper into the question of what exactly is the merger. 
First, the inspiral- merger-ringdown classification is based 
on the validity of the PN expansion and that of BHP. 
The inspiral, in other words, is just the set of times for 
which the PN expansion holds on the whole spatial do- 
main (to a given level of accuracy) . Correspondingly, the 
ringdown is the times for which BHP works everywhere 
throughout space. Merger, in this picture, is just the gap 
between those times during which PN and BHP theories 
give accurate results. 

In this paper, we propose that we can push each ap- 
proximation technique beyond its current range of use, as 
long as we do not apply it to all of space at a given time. 
We observe that at any time, there is a region outside 



a certain radius from the center of mass in which BHP 
applies. While this seemingly runs contrary to the com- 
mon notion that PN theory is the natural description of 
the weak-field region of a black-hole-binary spacetime, a 
black-hole metric in the limit of radii much larger than 
the mass and binary separation is identical to that of 
PN in the same region. If the PN expansion applies 
to the remaining portion of the spacetime (within the 
region where BHP holds), then BHP and PN theories 
could cover the entire physics of black-hole-binary coa- 
lescence. While it is somewhat unreasonable to suppose 
that PN theory truly applies to the strong-field region 
of a binary-black-hole spacetime, revisiting Price's treat- 
ment of non-spherical stellar collapse suggests that this 
may not be essential. 

In Price's 1972 paper [2l[, he addresses, among other 
issues, why aspherical portions of stellar collapse quickly 
disappear when, in fact, one could plausibly argue the 
contrary. Namely, if any non-spherical perturbations 
asymptote to the horizon (from the perspective of an 
observer far away), they would remain visible to this 
observer indefinitely. Price realizes, however, that the 
exterior of a collapsing star is just the Schwarzschild so- 
lution (due to Birkhoff's theorem, up to small perturba- 
tions), and perturbations on the Schwarzschild spacetime 
evolve via a radial wave equation in an effective poten- 
tial. Moreover, he notes that when the surface of the 
star passes through the curvature effective potential of 
the Schwarzschild spacetime, the gravitational perturba- 
tions induced by the star redshift. Finally, because the 
effective potential refiects low-frequency perturbations, 
the spacetime distortions produced by the stellar interior 
become less important, and observers outside the star ul- 
timately see it settle into a spherical black hole in a finite 
amount of time. Most importantly, this argument does 
not depend strongly upon the physics of the stellar inte- 
rior; as long as there is gravitational collapse to a black 
hole, this idea holds. 

In this paper, we adopt this idea, but we replace the 
stellar physics of the interior with a PN, black-hole- 
binary spacetime (see Fig. [T]). While in Price's case, the 
division of spacetime into two regions comes naturally 
from tracking the regions of space containing the star 
and vacuum, in our case the split is somewhat more ar- 
bitrary; one simply needs to find a region in which both 
PN and BHP theories hold, to some level of accuracy. 
How we choose the boundary between the two regions 
and the quantities that we evolve are topics that will be 
discussed in greater detail in Sec. HT] 

To test the above idea in this paper, we study a head-on 
collision of two black holes with transverse, anti-aligned 
spins and compare the waveforms and energy-momentum 
flux obtained from our approximation method with the 
equivalent quantities from full numerical simulations. 
Specifically, we organize this paper as follows. In Sec. [Hi 
we give a more detailed motivation for our model, and 
we then present the mathematical procedure we use in 
our method, for an equal-mass, head-on collision of black 
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FIG. 1: (Color online) This figure depicts a spacetime dia- 
gram of a black-hole collision, modeled after Price's descrip- 
tion of stellar collapse. We choose the trajectory of the two 
holes as a way to separate the spacetime into an interior 
and an exterior region. The exterior region is a perturbed, 
black-hole spacetime, whereas the interior is that of a post- 
Newtonian (PN) black-hole-binary system [shaded in yellow 
(light gray)]. The red (dark gray) region of the diagram shows 
the place at which the effective potential of the black hole is 
significantly greater than zero. This formalism allows us to 
divide the waveform into three sections: inspiral (or infall), 
which extends from the beginning of the binary's evolution 
until when the I = 2 effective potential of the exterior BHP 
spacetime starts to be exposed; merger, which extends from 
the end of inspiral to when the majority of the exterior poten- 
tial is revealed; and ringdown, which represents the remain- 
der. We overlay the even-parity, {I — 2,m = 0) mode of the 
waveform. 



holes. In Sec. IIIH we present an explicit calculation for 
the head-on collision of two black holes with transverse, 
anti-aligned spins, and we compare waveforms, radiated 
energy and radiated linear momentum, from our model 
with the equivalent quantities from a full numerical- 
relativity simulation. In Sec. IIVI we discuss how our 
method can help interpret the waveform during merger, 
and finally, in Sec.|Vl we conclude. We will use geometri- 
cal units {G = c = 1) and Einstein summation convention 
throughout this paper. 



II. A DETAILED DESCRIPTION OF THE 
METHOD 

A. Further Motivation 

Before going into the details of our procedure, it is 
worth spending some time discussing why our specific 
implementation of PN and BHP theories will help avoid 
some of the difficulties that arose in other methods in the 
introduction and noting the limitations and assumptions 
that underlie our approach. 

It is certainly hard to argue that existing orders of 
PN (up to v^ in the metric, for near-zone dynamics 12211') 
and BHP (up to second order for Schwarzschild, see [2^ 
for a gauge-invariant formulation) theories are accurate 
in the whole space, simultaneously. Nevertheless, it is 
plausible to argue that these approximation techniques 
cover different spatial regions at different times in a way 
such that each theory is either valid to a reasonable level 
of accuracy or occupying a portion of spacetime that will 
not influence physical observables where it fails. Using 
an approach of this type, we aim to get the most out of 
the approximation methods. 

Specifically, we find that the following procedure gives 
good agreement with the waveform of a numerical- 
relativity simulation presented in Sec. IIIII First, we have 
the reduced mass of the binary follow a plunging geodesic 
in the Schwarzschild spacetime. Then, we divide this tra- 
jectory in half to make a coordinate radius (and thus a 
coordinate sphere) that passes through the centers of the 
black holes. The set of all the coordinate spheres de- 
fines a time-like surface in spacetime. Finally, we apply 
PN theory within this time-like surface and BHP on the 
exterior. The two theories must agree on this time-like 
surface, which we will subsequently call the shell. 

Matching PN and BHP theories on this shell has cer- 
tain advantages. Because BHP theory relies upon a mul- 
tipole expansion, this makes it necessary to treat the PN 
interior in terms of multipoles of the potentials. For one, 
this is useful, because physical observables like the radi- 
ated energy and momentum very often do not need many 
multipoles to find accurate results. (In fact, in our exam- 
ple in Sec. Illli we see that the quadrupole perturbations 
alone suffice.) Second, a multipole expansion may also 
be helpful for the convergence of the approach. For two 
point particles, for example, each multipole component 
of the Newtonian potential U^ at the location of the 

particles satisfies Uj^ ^ M/R, where M is the mass of 
the binary and R is the distance from the center of mass. 
This is small for much of the infall, when R ^ M, and 
even when the binary reaches what will be the peak of the 
effective potential of the merged black hole, Uj^ ^ 1/3. 
Because the effective potential of the final black hole 
tends to mask perturbations within (as they are red- 
shifted and cannot escape), our hope (supported by the 
example in Sec. IIIip is that PN theory is still reason- 
ably accurate at the peak of the potential. Then, in our 



approach, the PN error will be hidden by this potential 
and, along with BHP in the exterior, it will suffice to 
explain the physics outside the black hole (in particular 
the gravitational waveform and the energy and momen- 
tum flux). Of course whether such a mechanism exists is 
not easily deduced analytically from first principles; only 
by comparing our results with those from full numerical 
relativity will we test this premise (which wc do in Sec. 
HIlBl . 

If this holds, our procedure has several advantages. For 
one, the matching works well at larger separations (which 
diminishes the influence of junk radiation). More impor- 
tantly, though, as the matching shell moves from large 
separations at early times to the vicinity of the horizon 
at late times, the spacetime smoothly transitions from 
predominantly PN to essentially BHP. One can see this 
most clearly in the way that the effective potential is re- 
vealed in Fig. [1] This leads to a waveform that smoothly 
transitions between an inspiral at early times to a ring- 
down at late times. Whether we correctly capture the 
merger phase is most easily confirmed by comparing with 
results from numerical relativity. Finally, our method can 
give a way to interpret the waveform during the differ- 
ent stages of a binary-black-hole coalescence by relating 
parts of the waveform with the retarded position of the 
matching shell at that point (for further detail, see Sec. 

There is also reason to suspect that the domain of BHP 
may extend to relatively early times in the merger, out- 
side the binary. A paper by Racine, Buonanno and Kid- 
der [2J] gives evidence in favor of this idea. In particular, 
they find, when studying the "super-kick" configuration 
(an equal-mass, circular binary of black holes with trans- 
verse, anti-parallel spins), that the higher-order spin- 
orbit contributions to linear-momentum flux dominate 
over the leading-order terms. These terms include both 
direct higher PN terms and tail terms, where the tail 
refers to gravitational radiation that is scattered off of 
spacetime curvature and propagates within the light cone 
(as opposed to the direct piece that propagates on the 
light cone). This suggests that even early on, the back- 
ground curvature plays an important role in generating 
the kick. Whether a Schwarzschild black hole properly 
represents this curvature is another idea that is difficult 
to argue for directly, but can be confirmed by comparison 
with the results of numerical relativity. 

In its current implementation, the feature of our ap- 
proach that is most arbitrary is setting the boundary 
between the PN interior and the BHP exterior. Never- 
theless, our procedure works for head-on collisions, and 
in future work, we will develop a more systematic way of 
treating the boundary. 



physics in the PN interior; (ii) we match the metric of 
the PN interior to the BHP metric through a boundary; 
(iii) we explicitly construct the boundary between the PN 
and BHP spacetimes; (iv) we evolve the metric pertur- 
bations in the exterior Schwarzschild spacetime; (v) we 
extract the waveforms and compute the radiated energy 
and momentum. Wc shall devote a subsection to each of 
these topics below. 

Before we do so, however, it will be helpful to clar- 
ify our notation regarding the different coordinates we 
use for the two metrics and the matching shell. In the 
PN coordinate system, we use Minkowski coordinates 
{T,X,Y,Z) and spherical-polar coordinates (i?, 0, $) 
within spatial slices. We will only consider linear per- 
turbations to Minkowski in the harmonic gauge. For our 
BHP, we employ Schwarzschild coordinates {t,r,9,ip)^ 
and similarly, wc only examine linear perturbations to 
the Schwarzschild spacetime. As we will show in Sec. 
Ill B 3| we can match these descriptions when R,r ^ M. 
This procedure is accurate up to terms of order (M/R)'^ 
in the monopole part of the metric and M/R in the 
higher-multipole portions, assuming that we relate the 
two coordinate systems by 



T = t, e = e, <^^ip,R = r-M. 



(1) 



The identification above allows us to label every point 
in spacetime by two sets of coordinates, (t, r, 0, ip) and 
(t, R,9,ip), where R = r — M. Because our program 
relies upon applying PN theory in an interior region and 
BHP on the exterior, it is therefore natural to talk about 
a coordinate shell at which we switch between PN and 
BHP descriptions of the space-time. On this shell, wc can 
use cither the Minkowski or Schwarzschild coordinates. 

In keeping with the notation above, we shall denote the 
separation of the binary by A{t) = 2i?(t) in PN coordi- 
nates and a(t) = 2r{t) in Schwarzschild coordinates. Fi- 
nally, we will denote the radial coordinate on the bound- 
ary by adding a subscript s to the PN or Schwarzschild 
coordinate radn [e.g., Rs{t) = A{t)/2 or rs{t) = a{t)/2]. 
We summarize the two coordinate systems and how they 
match, in Table Hand Fig. [2 



^.'' A(t)^ 2RJt) = 2(rJt)-M) ^ a(t) - 2l\ir 




R-^ PN 

Matching Shell BHP r=R+M 

FIG. 2: This figure shows, at a given moment in time, the 
Schwarzschild and PN radial coordinates, the binary separa- 
tion, and the position of the shell where we match the two 
theories. 



B. Details of the Implementation 

The procedure that we follow can be broken down, 
more or less, into five steps: (i) we describe relevant 



While we introduce a new method of matching PN and 
BHP theories, the idea of combining PN and BHP ap- 
proximations is not new. In fact, it is at the core of the 
Effective-One-Body formalism of Buonanno and Damour 



PN Spacetime 



Matching Shell 



Perturbed Schwarzschild Spacetime 



Coordinates (t,R,6,ip) {t, Rs(t), 6, if) or (t,rs{t), 6, ip) {t,r,6,if), r = R + M 
Binary Separation A(t) A(t) or a{t) a{t) 
Matching Radius R{t) = A{t)/2 Rsjt) = a{t)/2 - M ot rs{t) = A(t)/2 + M r{t) = a{t)/2 

TABLE I: This table summarizes our notation for the coordinates, the binary separation, and the matching radius that we use 
for the PN spacetime, the BHP spacetime, and the matching region between the two. 



[20|. In the EOB description, however, they match the 
point-particle Hamiltonians of PN and BHP theories, 
rather than joining the spacetime geometry. It would 
be interesting, as a future study, to see whether one can 
combine our procedure with that of the EOB to produce 
a geometrical EOB approach. Le Tiec and Blanchet [la| , 
on the other hand performed a more accurate matching 
between PN and BHP in their close-limit calculation with 
2PN initial data, but their matching only takes place on 
a single spatial slice of initial data. It would also be in- 
teresting to extend their higher-order approach to our 
procedure as well. 



1. The PN Interior Solution 

For our method, we will need to describe the metric of 
the PN spacetime in the interior, which wc do at leading 
Newtonian order: 



ds 



PN = -{l-2UN)dt 



+ {l + 2UN)6^, 



8w^dtdX' 
dX'dX^ . 



(2) 



Here Un is the Newtonian potential and Wi is the grav- 
itomagnetic potential, and the index i runs over X, Y, 
and Z. Our notation follows [25| [the above takes the 
results of Eq. (2.1) of that paper]. We then expand the 
Newtonian potential, Un-, into multipoles, keeping only 
the lowest multipoles necessary to complete the calcula- 
tion. In this paper the monopole and quadrupole pieces 
suffice (the dipole piece can always be made to vanish). 



Un^U^n''+U^n^ 



(3) 



potential, w, in spherical polar coordinates, we will be 
able to remove the radial component through a gauge 
transformation. Wc will, therefore, consider just the 9 
and if components of w, and when writing it in index 
notation, we will denote them with Latin letters from 
the beginning of the alphabet, (e.g. a,b = 9,ip). We can 



then expand the components Wa 
spherical harmonics. 



(2) 



in terms of two vector 



w, 



(2) 



w'^:T'^aY2,rn{0,V) 



'^bYo. 



(6) 



Here Vq is the covariant derivative on a 2-sphere, and 
Eq '' is the Levi-Civita tensor (with nonzero components 
ee^ = 1/ sm.0 and f-^p^ = — sin0). A convenient abbre- 
viation for the two spherical harmonics above is 






2.7n -M^2.m 



(7) 



The two terms arc denoted by (e) and (o) as a short 
hand for even and odd, because they transform as (—1)' 
and (—1)'"*"^ under parity transformations, respectively. 
The odd- and even-parity vector spherical harmonics are 
given explicitly by 



and 



A.O — 






1 dY^"^' 
sin 6 dip 



X^" =sin^ 



Qyl'. 



86 



dY 



Im 



89 



V 



dY 



Im 



dip 



(8) 



(9) 



respectively. These are the only parts of the PN metric 
that will be necessary for our approach. 



The quadrupole piece can be expressed as a term with- 
out angular dependence times a spherical harmonic, as is 
done below, 



U 



(2) 

N 



U^N"'Y2,miO,ip). 



(4) 



Wc shall follow the same procedure with the gravitomag- 
netic potential, although here we will, temporarily, keep 
the dipole term. 



w. 



(1) 



w. 



(2) 



(5) 



Wc will be able to remove the dipole term through gauge 
transformations, but this discussion is much simpler on a 
casc-by-casc basis. When we write the gravitomagnctic 



2. Matching to Perturbed Schwarzschild 

Wc then note that the Schwarzschild metric takes the 
form 



ds' 



-ll-'-^)dt'+(l-'-^] dr'+r'd^n, 



(10) 
where the last piece is the metric of a 2-sphere. We will 
use r without any subscript to denote the Schwarzschild 
radial coordinate. When M <C "r the Schwarzschild met- 
ric takes the form 



ds' 



- I 1 ]dt' 

r 



2M 
l + ±!i \dr'+r'd'n. (11) 



By making the coordinate transformation R = r — M , 
and identifying M/R with the monopole piece of the 

Newtonian potential, U}^ , then one can find that the 
Schwarzschild metric takes the form of the Newtonian 
metric in spherical coordinates, 



ds' 



r(0) 



(0)^ 



x{dR^ + R^d^n). 



(12) 



This similarity between the PN and Schwarzschild met- 
rics suggests a way to match the two at the boundary. We 
will assume that the monopole piece of the Newtonian 
potential becomes the M/r term in the Schwarzschild 
metric. For the remaining pieces of the Newtonian niet- 
ric (namely Uf^ and w,- ) we will translate them di- 
rectly into the Schwarzschild metric after performing 
any needed gauge transformations to make such a direct 
match reasonable. 

The original works on perturbations of the 
Schwarzschild spacetime are those of Regge and 
Wheeler [2a] for the odd-parity perturbations and Zerilli 
[27| for the even-parity perturbations. Moncrief [23| 
then used a variational principle to show that one 
can derive quantities from the metric perturbations of 
Regge, Wheeler and Zerilli that satisfy a well-posed, 
initial-value problem in any coordinates that deviate 
from Schwarzschild at linear order in perturbation 
theory. These quantities are related to the gravitational 
waves at infinity, and they satisfy a one-dimensional 
wave equation in a potential. We follow Moncrief's 
approach in computing these so-called gauge-invariant 
metric perturbation functions, but for our notation, we 
use that of a recent review article by Ruiz et al. |29j . 

Both the even-parity [transform as (—1)' under parity] 
and odd-parity [transform as (—1)'"''^ under parity] per- 
turbations are not very difhcult to find. By writing the 
PN metric in spherical-polar coordinates. 



r(0) 



r(2 



(0) 



r(2) 



ds' = -{l-2U'j:;'^2U'^')dt' + il + 2U'r:;'+2Ul,^') 
X [dR^ + R^d^n) - Sw^^^dtdx'' , (13) 

where dx'' = dO, dtp, one can see that the even-parity 
perturbations are nearly diagonal in the metric. In fact, 
at leading Newtonian order, it is exactly diagonal, be- 
cause the non-diagonal term coming from w/^ arises at 
a higher PN order. We will show this explicitly in Sec. 
mil For this reason, we only consider the diagonal metric 
components in the discussion below. 

The even-parity metric perturbations in Schwarzschild 
are often denoted 

[IT-tt J(o) - ^tt r , (ft^V j(o) - ^rr ^ ' (l^) 

\"'9B )(o) ^r K Y , {h^^ )(o) = r sm fA ■ y 

a specialization of Eqs. (57)-(59) of Ruiz et al. Thus, by 
matching the two metrics, one can see that 



The odd-parity term is somewhat simpler, because there 
is only one metric perturbation in Schwarzschild to match 
at leading order, 

{hts )(o) = h^ Xg , (V^ )(o) = V X^ , (16) 

Eq. (61) of Ruiz et al. From this, one can find that 



2,m 



-Aw 



2,m 

(o) ■ 



(17) 



The matching procedure thus gives a way to produce 
perturbations in the Schwarzschild spacetime. 



3. The Boundary Shell 

We now must find a boundary region where one can 
match a PN metric expanded in multipoles with a per- 
turbed Schwarzschild metric. For head-on collisions, we 
find that the boundary can be a spherical shell whose 
radius varies in time as the binary evolves. We can mo- 
tivate where this boundary should be with a few simple 
arguments, but the true test of the matching idea will 
come from comparisons with exact waveforms from nu- 
merical relativity. 

We know that at early times and for larger separations 
of the black holes, the PN spacetime is valid around the 
two holes; thus it is not unreasonable to suppose that 
the shell should have a radius equal to half the binary 
separation. Later in the evolution, BHP will be valid ev- 
erywhere, so the shell should asymptote to the horizon of 
the merged hole (as seen by outside observers). The tra- 
jectory of the shell should be smooth throughout the en- 
tire process, as well. Finally, the boundary should mimic 
the reduced-mass motion of the system, which physically 
generates the gravitational waves. 

A simple way to achieve this quantitatively is instead 
of having the motion of the reduced mass follow the 
PN equations of motion, we impose that it undergoes 
plunging geodesic motion in the Schwarzschild space- 
time. Given that the exterior spacetime is a perturbed 
Schwarzschild, and that we are matching the two approx- 
imations on a shell that passes through the centers of the 
two black holes, it is just as reasonable to use a trajectory 
in the Schwarzschild spacetime. Moreover, at large sepa- 
rations, the motion of the reduced mass of the system in 
both Schwarzschild and PN are quite similar; we primar- 
ily choose the geodesic in Schwarzschild for its behavior 
at late times. For completeness, we write down the dif- 
ferential equation we use to find the motion of a radial 
geodesic in Schwarzschild. Since we think of the black 
holes as point particles residing in the PN coordinate 
system, we write the evolution of the binary separation 
A(i), measured in the PN coordinates. 
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where i? is a positive constant {B = 1 ~ 2M/A{0) for a 
head-on plunge from rest). This expression can be found 
in many sources; see, for example [30| . The coordinate 
radius of the shell in the PN spacetime is just half the 
distance A{t), 



Rs{t)^lA{t). 



(19) 



Since the PN and Schwarzschild radii are related by 
R = r — M, one can find that in the Schwarzschild coor- 
dinates, the position of the shell is given by: 

r,it)^^Ait)+M. (20) 

In the Schwarzschild coordinates, since A{t) goes to 
2M at late times, the radius of the shell asymptotes 
to the horizon. More specifically, let us first define 
Srs{t) = rs{t) — 2M and then note that at late times, 
the trajectory approaches the speed of light as it falls 
toward the horizon. In terms of the tortoise coordi- 
nate, r, = r 4- 2Mlog[r/{2M) — 1], this occurs when 
V ~ t + r^ ~ const.. Writing this with respect to the 
variable 6rs{t), we find 
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We neglected the constant value of v because doing so 
has no affect on finding the scaling of Vg {t) at late times 
t. The equation above has a solution in terms of the 
Lambert W function, W{x), given by 
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Because we are interested in the behavior at large f, 
g-*/(2M) jg small, and we can use the leading-order 
term in the Taylor series for the Lambert W function, 
W{x) - X + 0{x^). We find that 

-t/(2A/) 
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Thus, matching the spacetimes at half the PN separa- 
tion of the binary and having the separation of the binary 
evolve via a Schwarzschild geodesic in the PN coordinate 
system makes the shell track the PN reduced-mass mo- 
tion at early times, but still head to the horizon at late 
times. We illustrate these different behaviors by plotting 
the full trajectory of the shell rs{t), the trajectory of a 
shell that follows a plunging Newtonian orbit (which we 
denote by 7's,(Af)(0 ^^^^ which represents the behavior of 
the shell at early times), and e"*/'^^-'^' (the late-time be- 
havior of the shell) in Fig. |31 We choose the initial values 
of the trajectory to conform with those of the numerical 
simulation with which we compare in Sec. lIIIBl The up- 
per and lower insets show how the shell's trajectory con- 
verges to the Newtonian behavior (at early times) and 
the expected exponential decay at late times. 

While choosing the trajectory of the shell might have 
a slightly ad hoc feel, in future work we will develop a 
framework that determines the shell motion consistently 
through radiation reaction. 
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FIG. 3: (Color online) This figure shows the trajectory of the 
boundary shell as the solid blue (black) curve labeled by Vs (t) . 
The other two curves show the early- and late-time behavior 
of the shell. The red (gray) dashed curve labeled by f^ (jv)(i) 
shows the trajectory of a shell that follows the Newtonian 
equations of motion for a plunge from rest. The green (light 
gray) dashed and dotted curve [denoted by e~'' "^^'] shows 
the exponential convergence to the horizon at the rate ex- 
pected in a Schwarzschild spacetime. The upper inset shows 
how the shell agrees with a Newtonian plunge from rest at 
early times, and the lower inset shows how the shell converges 
exponentially to the horizon at the expected rate. 



4- Black-Hole Perturbations 

One can then take the odd- and even-parity metric 
perturbations from the second subsection and transform 
them into two quantities, the Regge- Wheeler and Zerilli 
functions, respectively, that each satisfy a simple wave 
equation. We first treat the even-parity perturbations. 
Eqs. (63)-(65) of Ruiz et al. show how to take metric 
perturbations and transform them into the even-parity, 
gauge-invariant Zerilli function. Substituting our Eq. 
p^ into those three of Ruiz, we find 
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where dr is just the radial derivative with respect to the 
Schwarzschild radial variable. The odd-parity perturba- 
tions come directly from applying our Eq. (|17p to Eq. 
(67) of Ruiz et al. This gives that 
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The odd- and even-parity perturbations then evolve 
according to the Regge- Wheeler and Zerilli equations re- 
spectively. 



(5?-a.l+v'(L)M)*(c"o) = 0' 



(26) 



where r^, = r + 2M\og[r/{2M) — 1] is the tortoise coor- 
dinate. The potentials can be expressed most concisely 
via the expression 






(27) 



where X ^ l{l + 1) and 

[/' (r) = l, UUr) 



A(A + 2)r2 + 3M{r - M) 
(Ar + 3A/)2 ■ 

(28) 
1. These expressions 



Here A = (^ - l){l + 2)/2 = A/2 
follow Eqs. (5.3)-(5.6) of [H]. 

In our procedure, we find it easiest to evolve the Regge- 
Wheeler-Zerilli equations using a characteristic method. 
To do so, we define m = i — r, and w = < + r* and see 
that the evolution equation becomes 
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We will now discuss how we evolve our Rcggc-Wheeler- 
Zerilli functions with the aid of Fig. [TJ 

We must provide data in two places, the matching shell 
(in Fig. [T] it is the lower-left curve labeled by the points 
AefgD) and the initial outgoing characteristic (the line 
labeled by AB on the lower right). Once wc do this, how- 
ever, we can determine the Regge-Whcelcr-Zcrilli func- 
tions within the quadrilateral (with the one curved side) 
ABCD. We already discussed how we determine the 
shell in Sec. IIIB3| and the data we provide along that 
curve are just the Regge- Wheeler [Eq. (^5)) ] or Zerilli [Eq. 
p4l)] functions restricted to that curve. The data we 
must provide along AB are less well determined. If our 
computational grid extended to past null infinity, then 
we could impose a no ingoing wave condition. At finite 
times, we can still impose this boundary condition, but it 
leads to a small spurious pulse of gravitational radiation 
at the beginning of our evolution. To limit the effects of 
this, we keep the shell at rest until the junk radiation dis- 
sipates, and then we begin our evolution. At this point, 
the data along the line AB more closely represent those 
of a binary about to begin its plunge. 

With these data, we can then evolve the Rcgge- 
Wheeler- Zerilli equations numerically, using the second- 
ordcr-accurate, characteristic method described by 
Gundlach et al. in [31|. The essence of this method 
is that one can use the data at a point plus those at 
one step ahead in u and v, respectively, to get the next 
point advanced by a step ahead in both u and v. Ex- 
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plicitly, if one defines ^w = '^/.{u + Au,v + Av), 
*M/ = *};")(u + Aw,u), *£ = *|;")(u,w + Aw), and 



one has that 
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where Tc is the value of r at the center of the discretizcd 
grid. Because our boundary data do not always lie on one 
of the grid points in the (w, z;)-plane, we must interpolate 
the bottom point '^,'™.{u,v') to fall at the same value 

of V as the next boundary point at \l//™^(M -I- Au,v). 
As long as we do this interpolation with a quadratic or a 
polynomial of higher degree, it does not seem to influence 
the second-order convergence of the method. Finally, we 
can extract the Regge- Wheeler- Zerilli functions from the 
line BC in Fig. [T] as they propagate toward future null 
infinity. 



5. Waveforms and Radiated Energy -Momentum 

As we mentioned at the end of the previous section, it is 
easy to find the Regge- Whceler-Zerilli functions from the 
exterior of our computational grid; this is useful, because 
these functions are directly related to the gravitational 
waveform h, asymptotically. For radii much larger than 
the reduced gravitational wavelength, r ^ AGw/(27r), 
one has that 
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where _2y;m is a spin-weighted spherical harmonic. The 
above comes from Eq. (84) of Ruiz et al., which also 
contains a discussion about the spin-weighted harmonics 
in an appendix. One can then substitute this into the 
usual expressions for the energy and momentum radiated 
by gravitational waves. 
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(where n^ is a unit vector, a dot represents a time deriva- 
tive, and dfl is the volume element on a 2-sphere). A 
lengthy, but straight-forward calculation done by Ruiz et 
al. shows that 
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[their Eq. (91)]. For the components of the momentum 
we are interested in (in the xy-plane) combining their 
Eqs. (86)-(88) and (93) and using their definition P^^ = 
Px + iPy gives 
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The bar denotes complex conjugation. The coefficients 
aim. and bi^m are given by their Eqs. (41) and (42), wliicli 
we reproduce here 



Yi.m are the usual scalar spherical harmonics. One can 
then see that the nonzero coefficients of the spherical 
harmonics are 
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With the framework now in place, we are prepared to 
make a comparison with numerical relativity. 



III. HEAD-ON COLLISION OF SPINNING 

BLACK HOLES WITH TRANSVERSE, 

ANTI-PARALLEL SPINS 

In this section, we discuss the specific example of a 
head-on collision of equal-mass black holes with trans- 
verse, anti-parallel spins. We will specialize the general 
framework presented in Sec. |IT]to the current configura- 
tion in the first subsection and then make the comparison 
with numerical relativity in the second. 
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After applying the transformation of the PN and 
Schwarzschild radial coordinates, R — r — M, (and simi- 
larly A[t) ^ a{t)- 2M) one finds that 
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One can then substitute this into Eq. ([M)) to find the 
Zerilli function, 
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We have only kept terms to linear order in M /r in this 
calculation, because we only use Newtonian physics to 
calculate the gravitational potential. At the boundary of 
rs{t) = a{t)/2, the perturbation is constant at leading 
order, and varying only at higher orders. 



A. The Hybrid Model for the Head-on Collision 

We will mimic the configuration used in the numerical- 
relativity simulation for ease of comparison. We thus 
choose our two black holes, labeled by A and B, to have 
masses AIa — Mb = M/2, to start with initial separa- 
tion Xa = A{0)/2 = -Xb [A{0) = 7.8M in the numer- 
ical simulations and Y^ = Yb — Za — Zb = 0] and 
to have their spins along ±Z axis, respectively (so that 
S^ = 0.5M^ and S§ = -0.5A/| and all other compo- 
nents of the spins are zero). Though they initially fall 
from rest, as in the numerical simulation, we will denote 
their speeds by Va and Vb- 
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2. Odd-Parity Perturbations 

The calculation with the gravitomagnetic potential 
is slightly more difficult, because it involves additional 
gauge transformations. The gravitomagnetic potential is 
given by 



1. Even-Parity Perturbations 
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As we argued in Sec. |TT1 the even-parity perturbation 
will only rely upon the Newtonian potential, which has 
the familiar form. 
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Here Ra and Rb denote the distance from the centers 
of black holes A and B in the PN coordinates. We then 
expand the Newtonian potential, Um, into multipoles, 
keeping only the monopole and quadrupole pieces (the 
dipole piece vanishes). 
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These results appear, for example, in Eq. (6. Id) of [25| . 
The new symbols N\ and Ng represent unit vectors 
pointing from the centers of the two black holes in the PN 
coordinates, and V\ and Vj^ are the velocities of the two 
black holes. Expanding the gravitomagnetic potential to 
leading order in A{t) and simplifying the trigonometric 
portions of the equations below, we see that 
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The variables S and V are just the magnitudes of the 
spin and velocity of each black hole, respectively. For 
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this equal-mass collision, the spins have the same mag- 
nitude, and the velocities of the holes do, as well. We 
must then convert the gravitomagnetic potential into 
spherical-polar coordinates. 
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There is a dipole term in the component wr of the grav- 
itomagnetic potential, and this term will not evolve ac- 
cording to the Regge-Wheeler-Zerilli equation. One can 
remove it via the small gauge transformation to the met- 
ric, 



Finally, we make the transformation to the Schwarzschild 
radial coordinate, R = r — M (and similarly for A{t) = 
a(t) ~ 2Af ), to find that 
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We can then find the Regge- Wheeler function from Eq. 
(I25l) which is 
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As before, we keep only terms linear in M/r. At the 
boundary, the odd-parity perturbation is 
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where the bar refers to a covariant derivative with respect 
to the background metric (in this case flat space). Recall 
that the metric components, hu^ are related to the grav- 
itomagnetic potential, Wi, by hu = — 4wi. If we make a 
gauge transformation where the only nonzero component 
of S,n is 
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this has several important effects. For one, it eliminates 
htr, and it introduces a term. 
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into hu- This term, however, is of IPN order, and, since 
we are considering only the leading Newtonian physics, 
we will drop it. Then, it turns the remaining pertur- 
bation into the sum of odd- and even-parity quadrupole 
perturbations. Letting b = 9,ip^ one has that 
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If one were to include the even-parity, vector-harmonic 
term in the Zerilli function, one would need to take 
its time derivative. This means it enters as a ncxt-to- 
leading-ordcr effect, and we can ignore that term in our 
leading-order treatment. Thus, the relevant perturbation 
of the gravitomagnetic potential is 
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3. Energy and Momentum Fluxes 

Finally, because we only have quadrupole perturba- 
tions, the expressions for the energy and momentum 
fluxes greatly simplify. The energy flux, for the I = 2 
modes (taking into account that the m = ±1 are equal 
and the m = ±2 modes are equal, as well) becomes 
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and the momentum flux is given by 
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We have also used the fact that vj/ 
this head-on collision. 
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B. Comparison w^ith Numerical Relativity 

In this section we compare the results of our head- 
on collision of spinning black holes (with transverse, 
anti-parallel spins) with the equivalent results from 
a numerical-relativity simulation (see the paper by 
Lovelace et al. jSy for a complete description of the sim- 
ulation). Although the paper by Lovelace et al. dealt 
mostly with using the Landau-Lifshitz pseudotensor to 
deflne a gauge-dependent 4-niomentum and an effective 
velocity to help develop intuitive understanding of black- 
hole collisions, they also investigated the gauge-invariant 
gravitational waveforms and radiated energy-momentum 
(calculated from the gravitational waves at large radii). 
We will not attempt to study any of these Landau- 
Lifshitz quantities in this work, and, instead, we will just 
look into the gauge-invariant radiated quantities. Specifi- 
cally, for our comparison, we focus on the waveforms (the 
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I = 2 modes of the gravitational waves) and the radiated 
energy and momentum. 

In Figs. U and [51 we compare, respectively, the even- 

(2 2) 

parity perturbation ^cvcn and the odd-parity perturba- 



,(2,1) 



- even 

tion ^qj^'' from our method with the equivalent quan- 
tities from the numerical simulation SI featured in 
Lovelace et al. (Since the I — 2, m — —2,-1,0 com- 
ponents are related by constant multiples to the above 
perturbations, we do not show them.) In these figures, 
we also include a vertical, dashed red line that indicates 
the retarded time at which the shell crossed the light-ring 
of the final black hole r = 3M, in the hybrid method. 
This is the peak of the effective potential, and due to 
the low-frequency opacity of this potential, much of the 
influence of the boundary data is hidden within the po- 
tential after this time (and the waveform is due mostly 
to the quasinormal modes of the final black hole). Be- 
fore this time the match is not exact (as a result of junk 
radiation in the numerical simulation and the difference 
between the time coordinates), but the Newtonian order 
perturbations do quite a good job of exciting quasinormal 
modes of a reasonable amplitude. 
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FIG. 4: (Color online) Tlie blue (dark gray) dashed curve is 
^even from our hybrid method, whereas the black solid curve 
is the same quantity in full numerical relativity. The red (light 
gray) dashed vertical line corresponds to the retarded time at 
which the shell in the hybrid method reaches the light ring of 
the final black hole, r = 3M. We shift the numerical-relativity 
waveform so that the peaks of the numerical and hybrid 9 
waveforms align. 
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The even- and odd-parity waveforms in the hybrid 

(2 2) (2 1) 

method are the pieces of ^cvon and '^^^^ restricted to 
the outer boundary of the characteristic grid, labeled by 
BC in Fig. [TJ We found these perturbations through 
the procedure described in Sec. IIIB| applied to the spe- 
cific binary parameters described in Sec. IIII Al For the 
numerical-relativity waveforms, we chose to present them 
in terms of the even- and odd-parity perturbation func- 

(2 2) (2 1) 

tions ^ovcn and *o(j,i ' ^^ well. To find these perturba- 
tion functions from the numerical simulations, we first 




Numerical Relativity 
Hybrid Method 



-60 



-40 



-20 



20 

t/M 



40 



60 



80 



FIG. 5: (Color online) The blue (dark gray) dashed curve is 

f 2 1) 

'^odd from our hybrid method, whereas the black solid curve 
is the same quantity in full numerical relativity. The red (light 
gray) dashed vertical line corresponds to the retarded time at 
which the shell in the hybrid method reaches the light ring of 
the final black hole, r — 3M. We shift the numerical-relativity 
waveform so that the peaks of the numerical and hybrid ^^cven 
waveforms align. 



twice integrated the Weyl scalar ^4 with respect to time 
to get the waveforms /i+ and hx (since —^4 = h+ — ihx , 
at large radii, where a dot denotes a time derivative). 
One can relate them to the gravitational waveforms h^ 
and /ix by Eq- (EB), at large r. In the case of the 
/ = 2 perturbations shown here, rhj^'" = vG'I'cvcrn and 

rh^' = \/6*odd ■ ^® compared the /i+ and h^ found 
directly from the numerical simulation through extrac- 
tion of the Regge- Wheeler and Zcrilli functions from met- 
ric coefficients in the numerical code |33|, and the two 
procedures gave essentially identical results. 

In order to make the comparison between our hy- 
brid method and the full numerical-relativity waveforms, 
we must shift the numerical waveforms by a constant. 
Specifically, we choose this constant so that the peaks of 



(2 21 
the exact and approximate waveforms of ^ovon match (at 

a time that we set to be i = 0). We add this constant shift 
in time, because there is no clear relationship between the 
coordinate time at which the waveform in our code begins 
and the same coordinate time in the numerical-relativity 
simulation. Trying to find a relationship between these 
times is complicated by the fact that the hybrid method 
evolves on a characteristic grid, whereas the numerical- 
relativity simulation solves an initial-value problem in 
a gauge that changes as the black holes move together. 
Nevertheless, because both the numerical and the hy- 
brid method use asymptotically flat coordinates, at large 
radii, the time coordinates move at the same rate. This, 
in turn, means that it is only necessary to shift the time 
coordinates rather than rescaling them. As an interesting 
aside, Owen (SJ] found that this agreement between the 
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time coordinates in the numerical simulation and pertur- 
bation theory appears to extend even into the near zone, 
when he observed that multipole moments of the hori- 
zon oscillate at the quasinormal mode frequencies of the 
black hole. See the end of Sec. Ill of that paper for a 
discussion of why that might be the case. 

We also compute the momentum flux, and we show the 
accumulation of the velocity of the final black hole in Fig. 
[6l for both our method and the full numerical-relativity 
simulation. For our hybrid method, we use just the I = 2 
modes of the waveform to compute the momentum flux, 
Py, our Eq. (pT|) . We then find the velocity of the final 
black hole as a function of time by computing 



1 



yvit) = -JJ I Py{t')dt\ 



(62) 



to 



where we introduce an extra minus sign to account for 
the fact that the black-hole's velocity is opposite that of 
the momentum carried by the gravitational waves. For 
the numerical waveform, we show the equivalent velocity 
computed from the full Weyl scalar, ^4. For the numer- 
ical simulations, one typically computes 



Py — lim (/) sin 9 sin t 



■^Adt' 



dVL, 



(63) 



where ^4 is the Weyl scalar extrapolated to infinity, and 
drj is the surface-area element on a unit sphere. This ex- 
pression appears in a variety of sources [see, for example, 
Eq. (29) of the paper by Ruiz et al.]. We then can com- 
pute the velocity of the final black hole in the numerical 
simulations through Eq. (|62|) . as we did for the hybrid 
method. Again, we perform the same time-shifting pro- 
cedure as we did with the waveforms. The kick we find is 
remarkably close; 22 km/s for the numerical simulation 
and 25 km/s for our hybrid method. 

The radiated energy does not agree quite as well due 
to the fact that the even-parity perturbation is some- 
what larger than the equivalent numerical quantity (and 
it is the dominant contribution to the energy flux). Nev- 
ertheless, the results agree within a factor of two; the 
numerical simulation shows that roughly 0.057% of the 
initial mass is radiated, whereas our hybrid method pre- 
dicts approximately 0.096% of the initial mass escapes in 
gravitational waves. This is an improvement over some 
of the earlier, first-order, close-limit calculations at larger 
separations (see, for example, Fig. 1 of [35|, where more 
than 1% of the initial mass is radiated for equivalent ini- 
tial separations of the black holes). 

We are not aware of any equivalent calculations of ra- 
diated energy for head-on collisions within the EOB for- 
malism. While there have been recent 3PN calculations 
of the energy flux for a head-on plunge of black holes by 
Mishra and Iyer [3g, in PN-based calculations such as 
the EOB approach, one must always stop the calculation 
at some radial separation of the binary. For example, in 
Fig. 2 of Mishra and Iyer's work, they stop the calcu- 
lation at a point at which only roughly a tenth of the 
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FIG. 6: (Color online) The blue (dark gray) dashed curve 
is the velocity of the final black hole as a function of time 
(inferred from the gravitational radiation) from our hybrid 
method, using only the I = 2 modes of the wave. The black 
solid curve is the equivalent quantity in numerical relativity, 
computed from the full Weyl scalar, ^4. The red (light gray) 
dashed vertical line corresponds to the retarded time at which 
the shell in the hybrid method reaches the light ring of the 
final black hole, r = 3Af . As before, the numerical-relativity 
velocity is shifted so that the peaks of the numerical and hy- 



brid ^even wavefomis align. 



total energy that will be radiated in the head-on colli- 
sion, has escaped. This poses a small problem for EOB 
approaches, because as described in the introduction to 
this paper. Sec. U one must choose a point at which to 
stop the EOB inspiral-and-plunge waveform and match it 
to a set of quasinormal modes to obtain a full waveform. 
For inspiraling black holes, there is a natural point to do 
this: when the frequency of the inspiral-plunge dynamics 
approaches the quasinormal mode frequencies of the final 
black hole that will be formed. For a head-on collision, 
however, there is no analogous frequency at which one 
can match. We will, therefore, reserve any comparisons 
between our method and that of EOB for future work, 
when we extend our method to inspiraling, black-hole 
binaries. 



IV. THE THREE STAGES OF BLACK-HOLE 
MERGERS 



In addition to producing reasonably accurate full wave- 
forms, our approach also provides a possible interpreta- 
tion of the infall, merger and ringdown stages of a binary- 
black-hole merger. As shown in Fig. [1] before the shell 
reaches point e and enters the strong-field region [the red 
(dark gray) area, in which the 1 = 2, even-parity, effec- 
tive potential exceeds 1/3 of its peak value], the majority 
of the retarded waveform propagates along the light cone 
(the so-called direct part). The direct part overwhelms 
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the waves that scatter off the background curvature, be- 
cause, far away from the source, the curvature is small. 
This part of the waveform corresponds to the inspiral or 
infall phase. 

In the strong-field region, however, there is strong cur- 
vature (the black- hole effective potential in our model). 
While some fraction of the waves will pass directly 
through, as the shell enters this region at point e, waves 
that scatter off the effective potential (and thereby prop- 
agate within the light cone) become more significant. 
These waves often are called tail waves. Although PN 
waveforms do include the tail part, the fact that the 
higher-order PN terms that contain the tail dominate 
over the lower-order terms [Ijl does not bode well for 
the ability of the PN series to easily capture this ef- 
fect. Nevertheless, we are able to associate this mixture 
of direct and tail portions in the waveform to merger. 
In our model, this stretch of the waveform is related to 
the retarded times when the shell is passing through the 
strong-field region of spacetime (points e, /, and g in the 
diagram) . 

Finally, after the shell passes through the potential, the 
details of the perturbation no longer become important, 
as was found by Price in his stellar collapse model. Be- 
cause waves do not efficiently propagate through the bar- 
rier, the gravitational waveform associated with points g 
through D should arise from before and while the shell 
passes through the effective potential (not after). This 
last piece is that of ringdown. 

There is one subtlety to note about our interpretation 
of ringdown that might arise if the final black hole is a 
Kerr black hole. Mino and Brink [37| and Zimmerman 
and Chen [S^ showed that mergers that lead to a Kerr 
black hole can emit waves at integer multiples of the hori- 
zon frequency that decay at a rate proportional to the 
surface gravity. These modes come from a calculation in 
the near-horizon limit of a Kerr black hole, and from the 
vantage point of observers far away, these waves would 
appear to be coining from the horizon. These modes 
have a sufficiently high frequency that they could pen- 
etrate the effective potential of a Kerr black hole and 
contribute to the ringdown portion of the gravitational 
wave. Nevertheless, if we expand our description of the 
ringdown phase to include these horizon modes, our in- 
terpretation holds more or less as described above. 



V. CONCLUSION 

In this paper, we show, by examining the head-on col- 
lisions of spinning black holes, that a combination of PN 
and BHP theories gives a gravitational waveform that 
matches well with that of full numerical-relativity sim- 



ulations. We were able to do this not by applying the 
approximation methods to distinct times in the evolu- 
tion of the system, but by choosing regions of space 
in which the two methods work and finding that the 
waveform from black-hole-binary collisions can be pro- 
tected from lack of convergence in these approximations. 
Specifically, our method lumps all monopole pieces of a 
PN black-hole-binary spacetime into the Schwarzschild 
metric and treats the higher multipoles as perturbations 
of that Schwarzschild that evolve via a wave equation. 
Moreover, since PN and BHP theories describe the wave- 
form, this suggests that much of the nonlinear dynamics 
appearing in the gravitational waves of a head-on black- 
hole-binary merger can be well approximated by linear 
perturbations of the Schwarzschild solution. 

Our approach certainly cannot replace full numerical 
simulations. For one, we must test its validity for differ- 
ent kinds of coalescence by comparison with fully nonlin- 
ear numerical results. Nevertheless, we are hopeful that 
our method maybe be useful for gaining further under- 
standing of the spacetime of black-hole-binary mergers 
and for producing templates of gravitational waveforms 
for data analysis. To move towards these goals, we would 
need to make several modifications to our method (whose 
implementation we leave for future work) . Most of these 
changes revolve around finding a way to treat inspirals of 
black-holes binaries within our method. The most neces- 
sary addition would be finding a way to consistently treat 
radiation reaction within the formalism. This feature is 
essential for capturing the correct inspiral and plunge dy- 
namics. Also important for describing realistic physics of 
the ringdown would be to analyze the problem in a Kerr 
background. Each of these problems requires significant 
work, so we leave them for future studies. 
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